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Abstract. In this paper we present a nonconforming finite clement method 
for solving fourth order curl equations in three dimensions arising from mag- 
nctohydrodynamics models. We show that the method has an optimal error 
estimate for a model problem involving both (Vx) 2 and (Vx) 4 operators. 
The element has a very small number of degrees of freedom and it imposes the 
inter-element continuity along the tangential direction which is appropriate 
for the approximation of magnetic fields. We also provide explicit formulae of 
basis functions for this element. 



1. Introduction 

The magnetohydrodynamics (MHD) equations describe macroscopic dynamics of 
electrically conducting fluid that moves in a magnetic field. MHD model is governed 
by Navier-Stokes equations coupled with Maxwell equations through Ohm's law and 
Lorentz force. As an example, a resistive MHD system is described by the following 
equations: 

' p(u t + u • Vu) + Vp = i(VxB)xB + /iAu, 

V • u = 0,° 

B t -Vx(uxB) = -i(Vx) 2 B-^Vx((VxB)xB) 
-^(Vx) 4 B, 

V • B = 0, 

where p is the mass density, u is the velocity, p is the pressure, B is the magnetic 
induction field, 77 is the resistivity, 772 is the hyper-resistivity, p is the magnetic 
permeability of free space, and p is the viscosity. The primary variables in MHD 
equations are fluid velocity u and magnetic field B. 

MHD model has widespread applications in thermonuclear fusion, magneto- 
spheric and solar physics, plasma physics, geophysics, and astrophysics. Mathe- 
matical modeling and numerical simulations of MHD have attracted much research 
effort in the past few decades. Various numerical algorithms have been used in 
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MHD simulations; examples include finite difference methods, finite volume meth- 
ods, finite element methods, and Fourier-based spectral and pseudo-spectral meth- 
ods [27]- In [T2J [T31 Q31 [TSl [2T] , two-dimensional, incompressible MHD problems are 
studied in terms of finite element approximations of the stream function-vorticity 
advection formulation. Since MHD flow often develop sharp interfaces, adaptive 
/i- refinement techniques have been applied in MHD simulations 16, 25, 33J. Finite 
element computations of MHD problems in three-dimensions have been reported in 

El H7J |23l 12111321. 

In the existing finite element discreitzations for the above MHD model, a stan- 
dard pair of stable or stabilized finite element spaces are often used to discretize 
the velocity and pressure variables in the fluid equations. For the magnetic field 
variable B, however, at least two approaches are possible when the fourth order 
term (Vx) 4 B is not presented in the model, namely when electron viscosity r/ 2 = 0. 
One approach is to use the standard edge element ([23]) and the other approach is 
to use the Lagrange element after replacing (Vx) 2 B by — AB ( El [231 132] ) ■ Both 
these approaches will become more difficult when the fourth order term (Vx) 4 B 
is presented. We may still replace (Vx) 4 B by a biharmonic operator A 2 B. But 
discretizing a biharmonic operator in three dimensions is challenging. It requires 
220 degrees of freedom per element if a conforming finite element is used. One 
possible way to reduce the number of degrees of freedom is to use nonconforming 
discretizations which allow weaker inter-clement smoothness constraints but still 
provid convergent approximations. Among the class of nonconforming finite ele- 
ments for fourth order problems, Morley-type elements are special in the sense that 
they provide approximations with polynomials of minimal degree 18, 2S]. In |29j . 
a systematic construction of Morley-type elements is provided for solving 2m-th 
order partial differential equations in R™. In particular, we may apply the element 
in |29) with n = 3 and m — 2 consisting of piecewise quadratic elements to our 
system of biharmonic equations. This amounts to 30 degrees of freedom on each 
clement. This element provides a reasonable discretization of MHD equations when 
it is appropriate to replace (curl) 4 operator by the biharmonic operator. This ap- 
proach, however, may lead to difficulty for certain boundary conditions in practical 
applications. Indeed the treatment of boundary conditions is also an issue for the 
second order problem if (Vx) 2 B is replaced by — AB [10] . 

One more natural approach is to discretize the fourth order curl operator by 
some generalized higher order edge elements. But such type of edge elements are 
not available in the literature. The construction of such type of edge element is the 
main goal of this paper. 

Another possible approach to deal with the fourth order term (Vx) 4 B is to use 
operator splitting technique. Namely, one can introduce an intermediate variable 
a = (Vx) 2 B and then reduce the original problem to a system of second order 
equations. However, it is known that for some problems, such a technique cannot 
be applied. For example, when modeling the bending of simply supported plate on 
non-convex polygonal domains, the original biharmonic problem is not equivalent 
to the lower order system of two Poisson equations [21 [22]. In view of this, we 
consider discretizing the fourth order problem directly. 

In this paper, we investigate MHD equations that contain both fourth-order term 
and second-order term. In the literature, the major tool used for performing MHD 
simulations involving a fourth order equation has been the pseudo-spectral method 
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PQ. By choosing an appropriate formulation, we are able to construct a finite 
element approximation for this problem. This is a nonconforming finite element 
that involves only a small number of degrees of freedom. 

The rest of this paper is organized as follows. In Section 2, we describe a sim- 
plified model problem and the corresponding variational formulation. In Section 
3, we construct basis functions and provide the convergence analysis. Finally, in 
Section 4, we give some concluding remarks. 



2. Model Problem 

In the following, we introduce model problems for the fourth-order magnetic in- 
duction equations described above. Assume that fl C M 3 is a bounded polyhedron. 
By considering a semi-discretization in time and then ignoring the nonlinear terms, 
we obtain the following equations: 



(2.1) 



a(Vx) 4 u + ^(Vx) 2 u + 7U = f, in Q, 
V • u = 0, in Q, 



where V • f = 0, and the parameters a, /3, 7 > 0. We consider homogeneous bound- 
ary conditions, 

(2.2) u X n = 0, V X u = 0, on dtt. 

The above choice of boundary conditions arise naturally in the variational formu- 
lation given below. On the other hand, in the numerical simulations of the problem 
with pseudo-spectral method, one often uses periodic boundary conditions, e.g., 

mm- 

It is worth pointing out that the parameter a is usually much smaller than either 
/? or 7. This fact imposes some difficulties in designing robust numerical methods, 
as have been studied in the context of biharmonic problems, e.g., [20| [30] . 

The above fourth-order curl equations also arise from an interior transmission 
problem in the study of inverse scattering problems for inhomogeneous medium, 
e.g., 0. 

In order to provide an appropriate framework for our analysis, we define the 
following function spaces: 

iJ(curl;0) = {v G (L 2 (n)) 3 | V x v e (£ 2 (ft)) 3 }, 
# (curl; Q) = {v e ff(curl; Q) | v x n = 0, on dfl}, 

V = {v e i? (curl; fi) | V x v e H^{n)}. 
V is a Hilbert space with scalar product and norm given by 

(u,v)y = (V(V x u), V(V x v)) + (V x u, V x v) + (u,v), 

IMv = \/ (U, U)y . 

The following lemma gives a sufficient condition for a piecewisely defined function 
to be an element in V. 



Lemma 2.1. // v is piecewise smooth, v x n and V x v are continuous across 
element interfaces, then veF. 
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Using the following identity: 

(Vx) 2 u = -Au + V(V • u) 
and V • u = 0, the first equation in ( |2.1[ ) can be rewritten in the following form: 

(2.3) -aV x A(V x u)+^(Vx) 2 u + 7U = f. 

Multiplying Equation ( |2.3[ ) by the test function v and using integration by parts, 
we obtain the following variational formulation: 

(2.4) Find u G V such that o(u, v) = (f, v), V v G V, 
where the bilinear form a(-, •) defined on V x V is given by 

a(u, v) = a(V(V x u), V(V x v)) + /3(V x u, V x v) + 7(u, v). 

The well-posedness of the above variational problem follows from the Lax-Milgram 
lemma. 

The next lemma indicates that the weak solution satisfies the divergence-free 
constraint. 



Lemma 2.2. Assume V • f = 0, and let u be the solution of problem {2.4-). Then 
V • u = 0. 

Proof. Choose test function v = Vip where ip £ Cq°(0), then 

(u,Vp) = (f.Vp), 

hence, V • u = V • f = 0. □ 
3. A Nonconforming Finite Element 

In this section, we construct a nonconforming finite element to solve the fourth- 
order equation. One of the advantages for using a nonconforming element is that the 
number of degrees of freedom is small compared to that for conforming elements. 
The following construction is based on Nedelec elements of the first family that 
consist of incomplete vector polynomials |19j . The advantage of using incomplete 
vector polynomial space is that it provides the same order of convergence in terms 
of energy norms as the one given by corresponding complete polynomial space. In 
the following, we define the degrees of freedom in a special way to ensure that the 
consistency error estimate holds. 

Definition 3.1. The finite element triple (K, Px,^>k) is defined by 

• K is a tetrahedron; 

• Vk = Rz{K) = Pi © {p € (P 2 ) 3 | P • x = 0}, where P 2 is the space of 
homogeneous multivariate polynomials of degree 2; 
Sj<- is the set of degrees of freedom, see Figure [T] 

edge degrees of freedom: 



(3.1) 



M e {\i) = {Jvi-Tqds\VqePi{e), VeCifi, 

where r is the unit tangential vector along the edge e, 
face degrees of freedom: 



(3.2) M f (u) = ^(V x u) x n • q dA | V q € (P (/)) 2 , V / C xj, 

where n is the unit normal vector to the face /, 



A NONCONFORMING FEM FOR FOURTH ORDER CURL EQUATIONS IN R 3 5 

Ex- = M e (u)UAf/(u). 




Figure 1. Degrees of freedom of the finite element 

In the above finite element triple, the space Vk is the same as the second order 
Nedelec element of the first family for i/(curl) problem. The difference is the 
definition of the second set of degrees of freedom. It is designed specifically to 
ensure consistency for the fourth-order problems. The total number of the degrees of 
freedom for this element is 20, which is the same as the dimension of the polynomial 
space R-2,{K). 

It should be pointed out that the scaling factor 1/|/| 2 in the definition of the sec- 
ond set of degrees of freedom is associated with the construction of basis functions 
to be given later. 

The next lemma given in j!9] describes a relation between edge integrals and 
face integrals which will be useful in the error analysis. 



Lemma 3.2. If u € R,2{K) is such that the edge degrees of freedom (3.1) vanish, 
then 

(3.3) y (V x u) • n cL4 = 0, V face f <Z K. 



Proof. Given u £ R 2 (K) satisfies 

J u ■ r q ds — 0, V edge e C K. 

By Stokes' Theorem, 



/ (V/ x u T ) ■ q dA — I (V / x q) ■ u T dA = / u • t q ds, 
Jf Jf ' Jdf 

where ut is the tangential part of u, and V/x and V/X are surface vector curl 
and surface scalar curl, respectively. Let q be a constant. Notice that 



V/ x u T = (V x u) • n, 
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we conclude, 



J (V x u) • n dA = 0. 



□ 



As a direct consequence of Lemma 3.2, if both the edge degrees of freedom (3.1) 
and face degrees of freedom (3.2) vanish, then 



J (V xu)dA = 0. 



'/ 

The polynomial space R 2 (K) has the following property [8] . 
Lemma 3.3. If u e R 2 {K) satisfies V x u = 0, i/ien 

u = Vp, withp G _P 2 - 



We recall that the finite element (K, Vk, ^k) is said to be unisolvent if a function 
in Vk can be uniquely determined by specifying values for degrees of freedom in 
Ex. 



Proposition 3.4. The finite element defined by Definition \3.1\ is unisolvent. 
Proof. It is sufficient to prove that, given u £ R 2 (K), 

M e (u) = M f (u) = 0, Vec]f,/Cif^u = 0. 



Obviously, V(V x u) is a constant vector. Then using (3.3) and integration by 
parts, we obtain 

V(V x u) = -j- / V(V x u)dx = -j— I (V x u)n T dA = 0. 

\ K \ JK \ K \ JdK 

This implies that 

V(V xu)=0^>Vxu = const. 



Using again (3.3), we have 

V x u = 0. 



By Lemma |3.3[ we have 

u = Vp, withp £ P 2 (K). 

Since M e (u) = 0, we have 

J^qds^QyqeP^e). 
This implies dp/dr — on each edge e. Hence, p is constant and u = 0. □ 

In the following, we construct the basis functions. The explicit form of these 
basis functions not only is useful for implementation, but also instrumental for the 
interpolation error estimate. 



A NONCONFORMING FEM FOR FOURTH ORDER CURL EQUATIONS IN R 3 



7 




°1 



Figure 2. A tetrahedron with vertices <Zj, Oj, cifc, ai. and q^- 
are two tangential vectors on the face Fi. 

3.1. Basis functions. The main idea of the construction is to consider linear com- 
binations of basis functions of a related Nedelec element. Let K be an arbitrary 
tetrahedron with four vertices a^, a,j, a k and a/, see Figure [2j The corresponding 
barycentric coordinates are given by Ai, Xj, Afc, and A;, respectively. 

On each of the four faces, say face I (with vertices a,-,aj-,afe), we choose the 
following two tangential direction vectors: 

q tj = d~a\ = 6\K\(VX l x VA fe ), 

q 4 fc = aka\ = 6|if |(VAj x VA;). 
The edge degrees of freedom on edge (with vertices a t and cij) are defined 
explicitly by: 

Mj'(u)= / u-rds, 

where r is the unit direction vector of edge e^, s is an arc length parameter. The 
face degrees of freedom are defined as: 

M hj (u) = / (V x u) x n ( • q fj dA, 

\Ji\ Jf l 

M Hk (u) = / (V x u) x m ■ q lk dA, 

\fi\ Jh 

where n; is the unit outward normal vector of the face /;. 

We recall that the basis functions of the second order Nedelec element of the 
first family in barycentric coordinates are (see, e.g., [5J, [26], [31]): 

(1) Two basis functions on each edge e^-: 

Ljj = A^VAj — AjVAi, 

— A^VAj -|- AjY/Ai. 

(2) Two basis functions on each face /;: 

Lyfc = Ai(AjVAfe — AfeVAj), 
Ljife = \j{\V\k — AfeVAi). 
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In the following, we list a few useful facts about the geometry of a tetrahedron. 

(1) The unit outward normal vector of face /; is given by 

nvA,ir 

(2) The two tangential vectors of face fi are given by qy and q,^. 

(3) Let hi be the height of the tetrahedron corresponding to the face //, then 

|VA,| = i 
hi 

(4) Let | if | be the volume of the tetrahedron K, then 

W = I « ,( ^ X * )I = (^W^- 

Next, we construct basis functions in baryecntric coordinates. They provide a 
set of dual basis functions with respect to the prescribed degrees of freedom. 

Step 1. Construct eight basis functions {4>Uj} corresponding to the face degrees 
of freedom such that 

(3.4) MW(^) = 0, 
and 

(3.5) M mnp ((j)iij) — S m iS ni Spj . 

We use the basis functions of the second order Nedelec element as building blocks 
as they automatically satisfy the first condition (3.4 1. Using the facts listed above, 
we find that the basis functions corresponding to the facial degrees of freedom on 
face fi are given by the following: 

<t>Uj = 3\K\(Lnj - L yfe ), 

4>Uk = 3|lf|(Ljifc — Lyfc). 
By direct calculation, we have 

(V x L Hj ) x VA, • (VA, x VA fc ) dA 

[2A ( (VAi x VXj) + X l (VXi x VAj) - Aj(VA z x VA*)] 
[VA ; (V ; • V z ) - VA fc (VA, • VAfc)] dA 

j A, ; dAj [(VA, x VAj) • VA fc ](VA ; • VA,) 

J A, dAj [(VAj x VA,) • VA fc ](VA ; • VA,) 

r Jl %\K\hi 



and 



j (V x L ljk ) x VA ; • (VA, x VA fe ) dA = ~\fi\^ 
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Hence, 

M Uj {4> Uj ) 



\fl\ 



(V x (j) U j) x ni ■ q zj dA 



If! 



^ / Vx[3|if|(L Ki -L yfe )] 
l\ Jfi 



IVA/ 



\K\VXi x VA fe )cL4 



18^% 

l/il a 
18\K\ 2 hi 

\fi\ 2 



/ V x (L, y - Ujk) x VAj ■ (VA ; x VA fe ) dA 
Jfi 



Similarly, M Hk ((f> H j) = 0, and M f[ , ((j> Hj ) = where I' ^ I. 

Step 2. Construct twelve basis functions {ip^j |1 < i < j < 4, t = 1, 2} corre- 
sponding to the edge degrees of freedom such that 

(3-6) M$(V$) = *M»n,ii, 



(3.7) 



= 0. 



Here, we use the edge basis functions of the second order Nedelec element as 
building blocks since they satisfy condition (3.6 1. Since V x = 0, L,-j automat- 
ically satisfy condition (3.7 1. 

For functions Ly, we need to subtract from them a linear combination of face 
basis functions so that (3.6) and (3.7) hold. This can be done because by construc- 
tion, our face basis functions have no edge moments. This strategy for constructing 
basis functions can be found in US] . 

Finally, we can write the basis functions of the new element as the following: 

(1) Two basis functions on each face I (1 < I < 4): 

tplij = 3\K\(Liij - Ly fe ), 

4>uk = 3|if |(Ljift - Ljjfe), 

where L /y = A^AjVAfc - AfeVAj). 

(2) Two basis functions on each edge (1 < i < j < 4): 

4 X) = L ^> 



1 mnp (-^ij ) 0mnp ■ 

3.2. Convergence analysis. Let l~h = {Ki}^t 1 be a triangulation of the domain 
Q. On this triangulation we introduce the finite element space Vh and define the 
discrete norm || • ||^ by 

, 1/2 

l|v|U = E (IMIo,* + H V x v IIo,a + ||V(V x v)|| 2 ^) ' 
-Ken 

Consider the following discrete bilinear form: 
a h (u h ,v h ) 



E 

Ken 



a(V(V x u/,),V(V x v /l )) L 2 (A - ) + /3(V x u /l7 V x v h ) L 2 (x) 

+ 7(ufc, v h ) L 2( K ). 
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It is straightforward to verify that the bilinear form satisfies 

Ofc(v,v) > \\v\\l, VvG VJ,, 

|a h (u,v)| < ||u|| fc ||v|| h> V ueV + V h ,v eV h . 
The nonconforming finite clement discretization of problem (3.10) is: 
Find Uh € Vh, such that for all Vh £ Vh, 

(3.8) a h (u h , v h ) = (f, v h ). 

The convergence of the above finite element approximation can be analyzed through 
the following second Strang lemma |3] . 

Lemma 3.5. 

n I, / • t I, I, |a/i(u,w h ) - (f,w h )| 

||u — Uhll/, < inf \\u-v h \\h+ sup n n , 

v hev h w h ev h \\™h\\h 

where the first term on the right-hand side is called the interpolation error and the 
second term is called the consistency error. 

In order to estimate the consistency error we first define an average operator Pt 
on a face / by 

PfVf = T-rr / W dA. 
I/I Jf 

Since for any Vh € Vh, the quantity Jj V x dA is continuous, we know that Pf 
is well-defined for V x Vh- The following two lemmas are standard results. 

Lemma 3.6. Given any face f C K and w g (H 1 ^)) 3 , 

J |w - P f w\ 2 dA < h K \vr\{, K - 



Lemma 3.7. 



/ |w| 2 dA</^Vl6> + Mw|! )K . 



Next, we estimate the interpolation error and consistency error separately. 

3.2.1. Interpolation error estimate. Let K and K'^ be the two tetrahedra sharing a 
common face /, tk be the local interpolation operator for the second order Nedelec 
element of the first family, namely, given u € V , define r^-u such that 



rfj-u ■ t ds — u ■ t ds, V edge e C K, 



and 



y (r K u xn)- g dA = J (u x n) • q dA, V q e (P (f)) 2 , V face f <Z K. 



if J f 
Define u/ e such that 

M e (uj) = Af e (r A -u) = M e (u), 
M/(u/) = [M/(rjfu) + M f (r K ,u)}/2. 
If / C an, we set M/(uj) = A//(r K u). 
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Lemma 3.8. Given u 6 V , let u/ be defined as above, then 

||u-uj||&</i(|u|a + |Vxu|a). 

Proof. Let r^u be the global interpolation operator defined by r^u|j<- = r^u. By 
triangle inequality, 

||u - U/||h < ||u - r h u\\ h + \\r h u - m\\h. 
By the interpolation error estimate of Nedelec element, we have 

||u — r^ulU £ &(Ma + |V x u| a ). 
Notice that on each tetrahedron K, 



r K u - Uj = ^2 ^2 M mnp(ri<U - Uj) 
fci< 



J mnp: 



where {4> m np} ar e basis functions on face /, and {M mnp (-)} are degrees of freedom 
on face /. Using Lemma 3.7 and ||g np ||L 2 (/) = 0(h 2 ), we get 

2\M mnp {r K u - u 7 )| = \M mnp (r K > f u) - M mnp (r K u)\ 



< 



1 

17F 
i 



l/l 2 

< ||qnp|U»(/) , 



J (V x rjf/ u - V x r K u) x n ■ q np dA 

J (V x (r K ' f U - u) x n • q„ p cL4 + -jp J (V x (r K u - u) x n • q np dA 



|V x (rjf/u- u) x n|| L 2 (/) + ||V x (r K u-u) x n|| L 2 (/) ) 



~ l/l 2 4 

< h- 2 {h- l ' 2 \\V x (r K , f u-u)\\ 0iKU K' +/i 1/2 |V x (rji-u-uJlLjajjc) 

< ft-V 2 |V X U|2,JCUJC'. 

Notice VAi = O^ 1 ), and H^mnpllQ ^ — )) Cauchy-Schwarz inequality, 
we have 

1/2 / v 1/2 



||r K u - U/Ho.k < ^ ^ |M mn p(rjifU - r f u) 

< /i 3 |V x u| 2 ,s(jf), 

where S(if) = U^^T^^/n^aiC'. 
Hence, 



v/CK 



(3.9) 



''7, 



u - u/||o,n = (X] H 7 'A'U - u 7 ||^) 1/2 < /i 3 |V x u| 2)! 



A' 



By inverse inequality, we have 

|| V x (r h u - uj)|| ,n < ^ 2 |V x u| 2 ,n, 

||V(V x (r h u - u/))|| 0l n < J»|V x u| 2 ,n. 
Combining these estimates, we get 

\\r h u- ui\\ h < h\V x u| 2 ,n, 
and the desired estimate follows. 



□ 
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Remark: We note that the error estimate (3.9) indicates that r^u and u/ are 
super-close. Such type of estimate can not usually be obtained by the standard 
scaling argument (using Bramble-Hilbert lemma) . In our proof, we made use of the 
detailed information of the basis functions constructed in the previous section. 

3.2.2. Consistency error estimate. Given a tetrahedron K, in addition to the local 
interpolation operator tk, we introduce another local interpolation operator tk 
corresponding to the first order Nedelec element of the second family, namely, f^-u S 
{Pi{K)f CR 2 (K), and 

J ((fxu) ■ t) q ds — j (u ■ r) q ds, V q G -Pi(e), V edge e C K. 

Consider two tetrahedra K and K'^ that share a common face /. Given v/, £ Vh, 
define vj{ = VhIk- By definition, 



?kVk = rjc' Vic , on face /. 



Hence, 



(3.10) V / ip- [{r K v K ) x n] dA = 0, V <p e #(curl; fi), 

K JdK 

where n is the unit outward normal vector of dK. 
Consider the decomposition (see [5]): 



where div wk = 0, v?k ■ n|aK = 0, and G PziK). The following Lemma 3.9 
be found in 

Lemma 3.9. 

\\r K WK - Wi<-||o,K ^ h\\V x va-IIo,^- 



As a consequence of Lemma |3.9[ we have the following estimate. 
Lemma 3.10. 

\\r K w K - v K \\o,K ^ ^||V x Vjf ||o,ic ■ 
Proof. Using the interpolation operators defined above, we obtain 
Hence, 

rjcvjf - v K = f K ^K - wjf. 



By Lemma 3.9 we obtain 

\\tkVk ~ Vk \\o,k ^ ^||V x v K \\ 0>K . 



a 



Now, we can show the following lemma, which is critical for the consistency error 
estimate. 
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Lemma 3.11. For tp £ H(curl; ft), 



| [ V' K xn)dA< fe(IM|o,n + IIV x ^\\ , n ) ( ]T ||V 

K JdK \ K 



1/2 



X VfcH? K 



Proof. By the interpolation error estimates of the Nedelec elements 

||V x (f K v K - vk)||o,k < h\\V X Vk||i,K: 
Lemma [3.10| and Equation (3.10), we have 

V / ip ■ (v K x n) dA = V / <p- [(f K v K - v K ) x n] dA 

K JdK w JdK 



Y / (V x (p) ■ (r K v K - v A ') dx - cp ■ [V x (f A v A ' - v A )] dx 

K Jk 

< X!(H V X VIIo.kH^Va- - VA-|| ,A- + |M|o,A'|l V x (rKVK ~ V K )\\q,k) 

/ \ V2 



K 



< 



ft(IMIo,n + ||Vx v ||o,n) X)H Vxv *llU 



□ 



Next, we show the consistency error estimate for the nonconforming finite cle- 
ment approximation defined above. 

Theorem 3.12. Assume that u *E V is sufficiently smooth and v/j G Vt, then 



Mu, v h ) - (f,v h )| < h(\\V x A(V x u)i| + |V(V x u) 



|(Vx) 2 uj| + j|Vxu||) fc||Vx V/l ||0 



1/2 



Proof. By applying integration by parts, we get 
(V(V x u),V(V x v h )) K 
= -(A(V x u), V x v h ) K + (V(V x u) • n, V x v h ) 9K 
= -(V x A(V x u), v h ) K + (A(V x u), w h x n) 0A - 
+ (V(V x u) -n,V x v fc ) 8Jf , 



and 
Hence, 



(V x u, V x v h ) K = ((Vx) 2 u,v /l ) A - (V X u,v h X n) 9K . 

«h(u, v h ) - (f,v ft ) 
= ^ [a(A(V x u), v ft x n) aA + a(V(V x u) • n, V x v ft ) aA 

-/3(V x u,v ft x n) aif ] 
= [(aA(V x u) - /3V x u, V/l x n) 0A ] 



- [«(V(V x u) n,V x vfcjax] 
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By Lemma |3.11[ we have 

[(«A(V x u) - f3V x u, v h x n) dK ] 



KeT h 

< h(\\A(V x u)|| ,n + ||V x A(V x u)|| ,n + ||V x u|| ,o + || ( V x ) 2 u|| ,n) 

\ 1/2 



\ K 



X Vfc||?jf 



By Lemma [376] and the inter-element continuity of V x v/,, we get 
[«(V(V x u) -n,V x v h ) BK ] 



Ken 
< a 



J2 ( V ( V x u)-n-P/(V(V x u) -n),V x v h - P/(V x v h ))j 

KeTh fCdK 

< h\v(y x u)| ll0 fi £ l vxv ^iul • 



The theorem follows by combining the above estimates of the two boundary inte- 
grals. □ 

Finally, we have the following convergence result. 

Theorem 3.13. Let u and u/j be the solutions of the problems (3.10) and (3.16) 
respectively, then 

||u - Uftlo,/, + I V x (u - u ft )||o,ft + ||V(V x (u - Uft))||o,h < /i|| u lkn 

when u G (iJ 4 (ft)) 3 . 

Proof. Using the second Strang lemma, 

||u - u h | ,h + I V x (u - u/OIo.h + ||V(V x (u - u h ))\\ Q M 

< inf (|u-w/4o,/> + ||Vx(u-w h )||o,fc + |V(Vx (u-wji))| ,/i) 

a/,(u,w/i) - (f,Wfc) 

+ sup — , 

w h eV h ,w h ^o || V x w h \\i,h 

and previous lemmas, the desired inequality follows. □ 
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